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ABSTRACT 

The lifetime of the structure in grand design spiral galaxies is observationally ill- 
determined, but is essentially set by how accurately the pattern's rotation can be 
characterized by a single angular pattern speed. This paper derives a generalized 
version of the Tremaine-Weinberg method for observationally determining pattern 
speeds, in which the pattern speed is allowed to vary arbitrarily with radius. The 
departures of the derived pattern speed from a constant then provides a simple metric 
of the lifetime of the spiral structure. Application of this method to CO observations 
of NGC 1068 reveal that the pattern speed of the spiral structure in this galaxy varies 
rapidly with radius, and that the lifetime of the spiral structure is correspondingly 
very short. If this result turns out to be common in grand-design spiral galaxies, then 
these features will have to be viewed as highly transient phenomena. 
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1 INTRODUCTION 

Spiral arms are some of the most strikingly beautiful fea- 
tures of galaxies, yet even now we lack a full understanding 
of their origins. It is well known that these features cannot 
simply be interpreted as long-lived aesthetically-organized 
"material sp irals" of particularly bright stars - as stressed by 
lOortl Jl962l) . the strong differential rotation in disk galaxies 
would rapidly wind up and destroy any such arrangement. 
It was this "winding dilemma" that motiyjU^edjLm_^J^yiij 
(1966) to build on the earlier intuition of iLindbladl Jl95lf) 
and develop a theory by which spirals could persist as long- 
lived density waves, thereby apparently solving the problem. 
However, ever since this theory was first espou sed, questions 
have remained. For example, iToomrel il969T) argued that 
such spiral density waves would tend to dissipate on rela- 
tively short timescales and must be repeatedl y re-ex cited 
via a swing-amplifier, while Sell wood fc Kahnl jl99ll) pre- 
sented the interesting alternative "groove mode" mechanism 
for their initial excitation, in which spiral stru cture is pro- 
duced by a pinching of the disk. lSellwoo d ( 2000) has also ar- 
gued that almost all numerical simulations strongly support 
the notion that spiral structure is transient; indeed recent at- 
tempts to impose any kind of long-lived spiral density wave 
on a numerical simulation only se em to work if the spiral 
struc ture is very tightly wound up llYano. Kan-Ya fc Goudal 
2003), which is not what one finds in many real grand-design 
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spirals. From an observational perspective, it has also of- 
ten been noted that the most dramatic spiral galaxies like 
M51 all seem to have interacting companions, which sug- 
gests that such interactions play a role in continuously ex- 
citing an evolving spiral structure . This anecdotal impres- 
sion w as borne out in the study bv lElmegreen fc Elmegreenl 
(1983), which showed that grand design spirals are signifi- 
cantly more common in group environments where such in- 
teractions will be more frequent. 

There therefore remains the most basic question about 
the spiral structure that defines the appearance of so many 
galaxies: how long does it persist, or, put more evocatively, 
how long could extragalactic travellers stay away and still 
recognize their own galaxy when they got home? This ques- 
tion is of more than purely navigational interest, as the 
longevity of spiral structure has fundamental implications 
for our understanding of galaxy formation and evolution. 
The interpretation of the Hubble sequence depends strongly 
on the duty cycle of the open, grand design spiral structure 
that defines the late end of the sequence, while the timescales 
for dynamical evolutionary processes su ch as the heating 
of th e stellar component by spiral arms dJenkins fc Binnevl 
1990) also clearly depend on how long these structures last. 

The degree to which a pattern persists is dictated pri- 
marily by whether it has a single "pattern speed," Q p , the 
angular rate at which the pattern rotates. One of the dif- 
ficulties in studying this quantity is that it has no simple 
relation to the more directly observable physical velocity of 
material at different radii in the galaxy. The situation is 
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further complicated by the fact that only one component 
of this velocit y is directly accessible throug h the Doppler 
shift. However, iTremaine fc Weinberel l)l984 T hereafter TW) 
elegantly demonstrated that one can invoke the continuity 
equation to derive tt p from the observable material kine- 
matics solely by making the assumptions that the kinematic 
tracer is conserved, that it orbits in a single plane, and that 
the pattern rotates rigidly, so that f2 p is a single well-defined 
constant. Specifically, one can write the continuity equation 
as 



_ E+ -E^ + -E, H =0, 



(1) 



where T,{x,y,t) is the surface density of material orbiting 
in a plane with the Cartesian coordinates chosen to be 
aligned with the galaxy's observed major and minor axes. 
TW showed how one can integrate this equation twice to 
get rid of the unobservable velocity component and elimi- 
nate the numerically-unstable spatial derivatives, leaving 



O i: 



Tixdx — 



T>Vydx. 



(2) 



The terms in both integrals are all measurable (v y is just the 
observable line-of-sight velocity corrected for the galaxy's in- 
clination). These integrals can be determined for a range of 
cuts at different distances y from the galaxy's major axis. In 
principle, each such cut provides an independent measure- 
ment of fl p , but a more robust approach to calculating the 
pattern speed involves making a plot of one of these inte- 
grals suitable normalized, {x) = J Y>xdx/ J Eda; against the 
other with the same normalization, {v y ) — J Y,v y dx/ J Tidx: 
such a plot should yield a straight line passing throu gh the 
origin, whose slope is D. p (Merrificl d fc Kiiiikenlll99fih . 

The TW technique has now been successfully ap- 
plied to a number of barred galaxies, using stars as the 
tracer population and optical absorption line spectroscopy 
to measure their kinematics, with the interesting result 
that these objects do have we ll-define d pattern speeds and 
are rapidly rotating (^ ^ugrri. Dgbattigt|^fc Cgrshril 2003; 
lOerssen. Kuiiken fc Merrifieldl 1200.4 ICorsini| l2004h . How- 
ever, analyses of spiral galaxies have generally proved less 
straightforward for several reasons. First, spiral structure is 
geometrically a good deal more complicated than a bar. This 
more complex structure means that some parts of the struc- 
ture may be at an orientation in which the TW method has 
no discriminatory power (essentially because (a;) and {v y ) 
both become very small). It is also the case that some parts 
of the complex spiral structure are clearly transient features 
that take no part in any global pattern rotation, but can 
nonetheless produce significant spurious signal in the TW 
analysis. Fortunately, the fact that we can choose the values 
of y to which the analysis is to be applied means that we can 
to a large extent avoid the regions compromised in this way. 
A second issue is that obscuration and star formation in spi- 
ral arms mean that application of the continuity equation 
to the observable stellar component of these systems may 
not be valid (although may perhaps be possible using near 
infrared emission) . The technique has therefore usually been 
applied to radi o observations of the gaseous components, ei- 
ther in HI (e.g. lWestDfahilll998h or CO fe.g. lRand fc Wallinl 
2004), to avoid problems of obscuration. Formally, conver- 
sion of gas into stars means that the continuity equation 



is not applicable to this tracer, either, but star formation 
does not have the disproportionate effect on gas that bright 
young stars have on the optical emission, so the analysis is 
not unduly compromised. 

The most interesting problem in the application of the 
TW method to spiral galaxies, however, is that even after 
the above issues have been dealt with the method still seems 
to break down because plots of {x) against {v y ) do not yield 
straight lines for any choice of major axis position angle 
iDebattistal l)2003h and lRand fc Wallinl ()2004h explore in de- 
tail how such plots can be sensitive to this parameter]. In 
some cases, the departures from a line can be attributed 
to the limited quality of the data, but in others clear sys- 
tematic departures from the linear fit are apparent. By far 
the simplest interpretation of this failure is that these spiral 
arms do not have a single pattern speed, so their structure 
must indeed be evolving with time. 

In principle, one can use this breakdown of the TW 
method to constrain the spatial variation in Sl p and hence 
the lifetime of the current spiral structure. However, the 
original TW method was not particularly well formulated 
for such an analysis because the integrals in equation 
combine data from a wide range of radii, so they produce 
complicated spatially-averaged estimates of any varying pat- 
tern speed , whic h have no simple physical interpretation. 
IWestpfahll (Il998f) was able to use these values to derive the 
spatial variation in Q p under the assumption that it only 
changes slowly with position, but such slow variations are 
by no means guaranteed. He also presented a generalization 
of the technique that resulted in a tractable integral equa- 
tion if £l p only varies with position along the major axis of 
the galaxy, x, but since this choice of Cartesian coordinate 
arises simply from the viewing angle and not any intrinsic 
property of the galaxy, the pattern speed will not in reality 
depend solely on this coordinate. 

This paper therefore revisits the problem and derives a 
generalization of the TW method in which the pattern speed 
is allowed to vary arbitrarily rapidly in the radial direction 
within the galaxy. This is not the only possible way that 
O p could vary spatially, but it is the simplest physically- 
motivated option. It explicitly allows for the possibility that 
a galaxy may contain a number of distinct features at dif- 
ferent radii, such as bars and spiral arms, each with their 
own pattern speeds, and it also permits one to calculate 
a simple metric of the lifetime of a galaxy's current pat- 
tern. The remainder of this paper is laid out as follows. Sec- 
tion [5] presents the mathematical generalization of the TW 
method to a pattern speed that varies with radius, and Sec- 
tion|^]gives an example of its application to the grand design 
spiral galaxy NGC 1068. Finally, conclusions are drawn in 
Section ^ 



2 THE GENERALIZED TW METHOD 

Even if the pattern speed is allowed to vary with radius 
within a galaxy, it is still the case that the pattern at a 
given radius simply rotates around in azimuth <j> with time, 
so one can replace the temporal derivative in equation Q 
with 



<9E 

dt 



-n P (r) 



<jy. 



(3) 
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as in the original TW analysis. Making this substitution in 
equation Q, and integrating over all values of x to eliminate 
the unobservable v x component of velocity, we find 



Q p (r)^- dx - 



d 

— T,v v dx — 0. 
ay 



(4) 



Still following TW, we now integrate over y to eliminate a 
noise-sensitive spatial derivative, yielding 



Q p (r) dxdy + 



T,v y dx — 0. 



(•>) 



Changing variables from Cartesians to polars in the double 
integral gives 



-arcsin(y / r) 



9E 

£l p (r)— — r dr d<j> + 



T,v y dx = 0, (6) 



/ r=y J (p— arcsin(y/r) ^ 

and carrying out the integration with respect to 4> then 
yields the final integral equation, 

poo 

{[V{x',y)-V(-x',yj\r}n p {r)dr= / Xv y dx, (7) 



smooth function that best matches the two sides of equa- 
tion 0. However, for this initial test we will avoid imposing 
any such prejudice upon the answer, as there is always a 
danger that a regularized solution owes its form more to 
the prior assumption about the smoothness of the function 
than the unpredictable astronomical truth. This freedom is 
particularly valuable in studies of pattern speed, where spa- 
tially distinct features can have completely different pattern 
speeds, leading to possible discontinuities in Q p (r). 

Once the variation in pattern speed with radius has 
been calculated from the solution of equation JHJ, it is simple 
to evaluate how the current spiral pattern will evolve into 
the future. By following this shearing, different metrics of 
the lifetime of the spiral structure can be constructed. One 
straightforward model-independent measure is provided by 
calculating the dispersion in the derived values of £l p (ri), 
<j(Q p ); a measure of the lifetime of the pattern is then just 
2n/a(il p ). Even simpler, one can just take the largest and 
smallest values within some range of radii, and calculate a 
characteristic wind-up time for this part of the pattern of 



Twind = 27r/(Sl™ ax — fi™ ln ) 



where x'(r,y) = -J V 2 — y 2 . As a check, it is straightforward 

to show that equation @ reduces to equation © in the case 3 APPLICATION TO NGC 1068 



These coordinates can be straightforwardly constructed 
from the observed coordinates, allowing for the galaxy's in- 
clination: x — x b s , y = y bs/cosi and v y — u b s /sini. 
Since the kernel term in curly brackets and the integral on 
the right hand side of equation J7J are both observationally- 
determined quantities, the problem is now reduced to 
a Volterra integral equation of the first kind for fi p (r). 
Such equations a re straightforward to solve numerically 
iPress et al.ll992l chapter 18): reduction of the integral to a 
discrete quadrature for different values of r — Ti and y — yj 
turns this problem into a matrix equation of the form 



E 

ri>Vj 



K{ri,yj)Q.p(ri) = f(yj). 



(8) 



The infinite upper limit on the integral complicates the anal- 
ysis a little, but as for the original TW method one finds 
that the decline in E with radius and the lack of strong spi- 
ral structure at large radii mean that the kernel decreases 
rapidly to zero, so that the integral has converged at an at- 
tainable radius. The lower limit on the integral means that 
the matrix K is triangular, and hence that equation J5J is 
readily soluble by back substitution to obtain Q p (ri). How- 
ever, the triangular form of the matrix does mean that a 
problem with the data at some radius will propagate in- 
ward to compromise all the derived pattern speeds at smaller 
radii, so a degree of care must be taken in evaluating the 
range of radii over which pattern speeds can be determined 
reliably. Note, though, that in principle one gets two in- 
dependent measurements of £l p (r), doubling the chances of 
finding a region free of such issues: here we have assumed 
that y > 0, but we could equally well derive the function 
from the data at y < 0. 

As formulated here, Q p can vary arbitrarily with radius. 
One could impose a degree of smoothness by regularizing the 
solution through an algorithm that penalizes rapid fluctua- 
tions in Q p (r), or even by finding a suitably-parameterized 



As an initial experiment, we have applied this algorithm 
to the BIMA Su rvey of Nearby Galaxies (BIMA SONG; 
iHelfer et al1l2003ft CO observations of the grand-design Sb 
galaxy NGC 1068, which offers a test case well suited to 
this analysis. Although optic ally classified as a normal spiral 
iSandaee fe Tammannlll987l) , albeit one with a Seyfert nu- 
cleus , the 2MASS near-infrared observations llJarrett et alJ 
2003) reveal a clear central bar, which provides a strong can- 
didate for the driver of the grand design spiral seen in the CO 
emission. iRand fc Wallinl ||2004) showed that NGC 1068's 
gas content is dominated by molecular material and that 
there is not much current star formation, so little of the 
molecular phase is being converted to or from other phases 
and the application of the continuity equation to this com- 
ponent alone is valid. 

A further concern in the application of the continuity 
equation to these data is how well the observed CO emis- 
sion traces the molecular gas surface density, and hence 
whether the equation can be applied directly to the data. 
In particular, even within the Milky Way it is known that 
the conversion between CO emission and molecular gas col- 
umn, c onventionally termed X, sh ows a metallicity depen- 
dence IBrand fc Wouterlootl ll995). while there is also evi- 
dence that X varies linearly with metallicit y from galaxy 
to galaxy ( Bos^h^^Lecrueux^^^Iayazzi 2002). However, as 
can be seen from iPilvughi^Vncn^z^^Contmil <|2004|) . there 
is essentially no metallicity gradient in NGC 1068 over the 
range of radii probed by the data analyzed here, so X should 
be close to constant. There is also the possibility that X 
could vary with other parameters. For example, enhanced 
macroscopic opacity effects, more intense radiation fields, 
or stronger deposition of CO on to dust grains could ren- 
der the values of X in spiral arms systematically different 
from those in lower density regions. However, recent stud- 
ies of molecular cloud populations in nearby galaxies with 
higher resolution and sensitivity than previously achieved 
are not finding any evidence for significant variations in 
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Figure 1. Plot of (s) versus (u,,) for the BIMA SONG CO ob- 
servations of NGC 1068. The numbers indicate the distance y t, a 
in arcseconds by which each cut is offset from the galaxy's major 
axis, and the dotted line joins adjacent cuts. 



X wi t h environment iRosolowskv et alj|2003t : iBolatto et alJ 
120031: iLerov et alJ 120051) . Physically, one might expect X 
to vary with temperat ure and density as X oc Tp~ ' 5 
jYoung fc Scovilldll99lft : the fact that there is so little ob- 
served variation then implies that in the astrophysical envi- 
ronment this combination of temperature and density varies 
very little in regions where there is a significant amount of 
molecular material. Interestingly, even if one does allow X 
to vary by a fairly generous factor of two between the arm 
and interarm regions, it makes very little difference to a TW- 
style analysis JZimmer. Rand fc McGrawl2004l) . Clearly, the 
matter has not yet been finally settled, and any conclusions 
will have to carry this caveat, but currently there is no rea- 
son to believe that variations in X should compromise this 
analysis. 

What makes NGC 1068 such an interesting target for 
this analysis is that the conventional TW analysis produce d 
an interesting result for this galaxy jRand fc Wallinll2004) . 
As Figure Q reiterates, although the plot of (x) versus {v y } 
roughly describes a straight line, it is apparent that there are 
significant departures from the expected relation. It is also 
clear that these departures are not the random noise that 
one would expect if they simply represented a limited signal- 
to-noise ratio - ordering the points by their distance y b B 
from the galaxy's major axis produces a systematic figure- 
of-eight in this plane. The deviant point at y b s = —30 arcsec 
arises due to a bright spot in the CO emission that this cut 
intersects, which is clearly not a part of the overall spiral 
structure, so, as discussed above, must be avoided in the 
TW analysis. With the exception of this point, however, 
NGC 1068 seems to show evidence for a systematic variation 
in pattern speed with position, warranting further investi- 
gation through the analysis developed in this paper. 

A preliminary indication of the variation in pattern 
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Figure 2. Plot of the weighted mean pattern speed as a function 
of distance from NGC 1068's major axis. Filled symbols are for 
positive y and open symbols show the corresponding values for 
negative y. The scale conversion from arcseconds to kiloparsecs 
has been made by adopting a distance to NGC 1068 of 14.4 Mpc 
feland-Hawthorn et al]|l997l) . 



speed with radius can be obtained by calculating (Q p ) = 
{v y )/{x) for cuts at different values of y bs- Figure |2] shows 
this quantity for both positive and negative values of y h B - 
Both sides of the galaxy show a significant decline in 
(flp) with radius, which, as discussed above, implies that 
NGC 1068 cannot have a constant pattern speed. In fact, 
since the value of (fi p ) at any given value of y b s is a weighted 
average of the pattern speed at all radii r ^ y, this plot will 
if anything tend to wash out the variation with radius, so we 
might expect the physical Q p (r) to be an even more steeply 
declining function. It is also notable that the two sides of 
the galaxy show a systematic difference in (fi p ), which pro- 
vides further evidence that the galaxy ca nnot be de s cribed 
by a global pattern speed - interestingly, IWestpfahll (119981) 
found a similar difference in his analysis of the spiral struc- 
ture in M81, suggesting that this phenomenon may be quite 
common. It should also be noted that the deviant point at 
2/obs = —30 arcsec does not show up dramatically in this 
plot, as both (a;) and (v y ) are discrepantly large, so their 
ratio comes out more-or-less in line with the other values. It 
is only by using the extra information in Figure that the 
departure from the grand design structure becomes appar- 
ent. 

Bearing these points in mind, it is still worth seeking to 
apply the generalized TW analysis to see how rapidly the 
pattern speed would have to vary with radius to be con- 
sistent with the observations. The fact that the two halves 
of the galaxy show somewhat different properties indicates 
that the evolution of the pattern must be somewhat more 
complex than this radial shearing, but it should nonetheless 
give a reasonable first approximation, and an estimate of the 
timescale involved for the pattern to evolve significantly. 
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Figure 3. Plot of the variation in pattern speed with radius for 
NGC 1068, as derived using the generalized Tremaine- Weinberg 
method. The panels show the results obtained using increasingly 
finely-binned discretizations of the integral equation J7J to form 
equation JHJ- 



To this end, we have solved the matrix equation (JSJ for 
the data at y > (thus avoiding the discrepant point at 
negative y) to estimate £2 p (r). To test the robustness of the 
result and estimate the associated errors, Figure|3]shows this 
analysis repeated using different radial binnings to discretize 
the integral equation (J7J into the matrix equation (|HJ . Even 
for these "well-behaved" data, the matrix inversion method 
fails inside ~ 20 arcsec. At these radii, the galaxy becomes 
close to axisymmetric, making the integrals in equation J7J 
small and meaning that noise starts to dominate. The fur- 
ther amplification of this noise by the inversion process then 
leads to wildly oscillating random values of Q p , a clear sign 
of numerical instability. However, at larger radii the differ- 
ent binnings produce very consistent rather smooth answers, 
indicating that the solution to equation J7J is numerically 
stable in this region even without regularization. The varia- 
tions between different binnings give a good indication of the 
uncertainty in £l p (r): although there are variations between 
the different binnings, the overall decline in Q p with radius 
is robustly reproduced. As Figure [5] already suggested, the 
pattern speed decreases strongly from ~ 150km/s/kpc at 
r ~ 1.5 kpc to ~ 50km/s/kpc at r ~ 2.5 kpc. This analy- 
sis fits surprisi ngly well with an entirel y independent study 
undertaken by ISchinnerer et alJ B200CI) . who looked at the 
likely locations for resonances in this galaxy and suggested 
that if an inner bar existed at r ~ 14 arcsec then it should 
have a pattern speed of ~ 140km/s/kpc, and that an outer 
oval-shape distortion which reaches out to r ~ 120 arcsec 
may have a pattern speed as low as ~ 20 km/s/kpc. The ad- 
vantage of the current study is that we are measuring pat- 
tern speeds directly without assuming anything about the 
resonant radii, such as presuming that the outer oval struc- 
ture in this system ends at its co-rotation radius. Figure [3] 



provides direct evidence that the grand-design spiral struc- 
ture between the inner bar and outer oval is not locked to 
either of their pattern speeds, but rather that the spiral arm 
pattern speed varies continuously with radius between these 
limiting values. This strong shearing implies a short lifetime 
for the present spiral structure between r ~ 20 arcsec and 
r ~ 40 arcsec of only r w i n( j ~ 10 s years, which is directly 
comparable to (but actually a little shorter than) the corre- 
sponding orbital period at these radii. This beautiful grand 
design spiral would appear to be only a transient feature on 
the face of NGC 1068. 



4 CONCLUSIONS 

This paper presents a generalization of the Tremaine- 
Weinberg method for calculating pattern speeds of galaxies 
which accommodates the possibility that the pattern speed 
may vary with radius. This generalization has potential ap- 
plications to a variety of different galactic systems. One 
could, for example, attempt to determine the distinct pat- 
tern speeds of multiple featu res within a single galaxy, such 
as ne sted bars-within-bars dCorsini. Debattista fc Aguerril 
120051) . 

Here, however, we have investigated the apparently sim- 
pler case of the grand design spiral structure in the CO emis- 
sion from NGC 1068. The conclusion of this analysis is that 
this spiral structure does not have anything close to a single 
pattern speed, and is winding up on a timescale marginally 
shorter than the orbital period, reintroducing the winding 
dilemma with a vengeance. In fact, we know that some of the 
other assumptions made in the analysis are violated at some 
level: the difference between the kinematic signal at positive 
and negative values of y implies that the pattern cannot 
be evolving with time in a manner entirely consistent with 
the simplest azimuthal winding up considered here. How- 
ever, this extra complexity simply means that the current 
design must be distorting even more profoundly and rapidly 
than we have found. In the only other comparable study 
to date. I Westpfahl jl998) tentatively reached a similar con- 
clusion regarding the grand design spiral structure in M81, 
based on both differences between the kinematics of the two 
sides of the galaxy and the inconsistency of the results from 
a TW-style analysis. If this result turns out to apply com- 
monly to other galaxies, then intergalactic travellers would 
be well advised not to use the morphology of spiral structure 
to identify their homes. 
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